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ABSTRACT 



The effect of the Monterey Submarine Canyon on seiching in Monterey 
Bay is not well known. Spectral analyses of simultaneous tidal records 
from the north-south extremities of the bay were performed for 23 
January and 20 April 1969 to investigate this effect. Both day's 
records had long-wave activity of which seiching was at least a con- 
tributing mechanism. Analyses of the computed spectra for the periods 
during the long-wave activity, and ten-hour periods both before and 
after, indicate that the seiching motion in Monterey Bay has similar 
amplitudes at the north-south extremities. 
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I. INTRODUCTION 



Monterey Bay, a large, semi -elliptical bay, is located approximate- 
ly sixty nautical miles south of San Francisco. It is bounded on the 
north by Point Santa Cruz and on the south by Point Pinos. The bay 
is bisected by the Monterey Submarine Canyon. The head of this deep 
canyon is located just off Moss Landing (Figure 1). 

Several studies have been conducted to investigate the long- 
period waves in the bay and its adjoining harbors. (For purposes of 
this research, a long-period wave is defined as any wave with a period 
greater than 1 min.) Hudson D949H investigated the surge in Monterey 
harbor for the U.S. Army Corps of Engineers. The study provided data 
for a proposed model of the harbor. He used six months of continuous 
data (October 1946 to April 1947) obtained by three automatic Stevens 
water-level recorders which were electrically synchronized. The 
recorders were located within the present boundaries of the harbor. 
Hudson's results for long-period waves are shown in Table. I. The 
cause of these long-period waves in the harbor was not known and 
Hudson did not attempt to define the surge mechanism. 

Wilson, et al D965II , also working under the auspices of the 
Corps of Engineers, made a field study for a proposed surge-action 
model of Monterey harbor. Wilson analyzed the long-period waves 
in the entire bay, using data from several sources. Wave recorders 
installed by Marine Advisers (MA) at Monterey and Santa Cruz, the 
north-south extremities of the bay, provided excellent data. Three 
sensors were installed within the harbor at Monterey and were in 
continuous operation from October 1963 to April 1964. Two were 
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arranged so that the tides and sea-swell were filtered out; a third re- 
corded the sea-swell approaching the harbor. Wilson found no correla- 
tion between the long waves in the harbor and the approaching sea-swell. 
This finding led him to conclude that the surge in the harbor is not a 
result of the incoming swell or surf beat. 

The MA sensor in Santa Cruz was inside the harbor, CGrauzinis, 19681]. 
Again, the short-period waves and swell were filtered out so that it 
functioned as a long-period wave recorder. This sensor was in con- 
tinuous operation from October 1963 to February 1964. A summary of the 
Monterey and Santa Cruz data is presented in Tables II and III. 

To obtain an independent evaluation of the oscillations in the bay, 
Wilson performed residuation analyses on several different records for 
various locations around the bay. Residuation analysis is accomplished 
by successively subtracting apparent periodicities from a wave record. 

The procedure is continued until a relatively smooth trace remains, 
and a sequence of apparent periods of oscillations is obtained. These 
results are presented in Table IV. Table V presents a synopsis of 
spectrum analyses for three days record of the Monterey MA sensors. 

In order to determine the cause of the surge, Wilson made an ex- 
tensive analysis of the two- and three-dimensional (i.e. two spatial 
and one time dimensions) oscillating characteristics of the bay using 
(1) approximate analytical solutions for a semi-enclosed basin, (2) 
numerical solutions for the modes of both two-dimensional and three 
dimensional oscillations, and (3) wave refraction diagram techniques. 

He concluded, based on the best fit of the observed sequences 
of periods in Tables IV and V with calculated modes of oscillation 
for the bay, that the Monterey Submarine Canyon causes the bay to 



13 



TABLE I 



Long-Period Waves in Monterey Harbor (after Hudson, 1949) 



PERIOD 


AVERAGE HEIGHT 


PERCENTAGE OF 


(min) 


(ft) 


Time PRESENT 


1-2 


0.4 


20 


2-4 


0.5 


30 


4-15 


not given 


15 



TABLE II 



Marine Adviser's Data for Monterey Harbor (after Wilson, 1965) 



SENSOR 


PERIOD 


AVERAGE HEIGHT 


PERCENT OF 




(min) 


(ft) 


TIME PRESENT 


1 


1.7-14 


0.1-2. 5 


0-50 


2 


1-14 


0. 1-3.0 


0-55 



TABLE III 

Marine Adviser's Data for Santa Cruz Harbor (after Wilson, 1965) 



PERIOD 


AVERAGE HEIGHT 


PERCENT OF 


(Min) 


(ft) 


TIME PRESENT 


1-2.3 




0-18 


2.3-4 


0. 2-2.0 


20-70 


10-14 




0-25 



14 
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Periods of Oscillation from Residuation Analysis (From Wilson, 1965) 
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Periods of Oscillation from Spectral Analysis (From Wilson, 1965) 



essentially act as independent north and south open-mouth oscillating 
basins with the boundaries lying along the natural mouth of the bay apd 
the center line of the canyon. That is, oscillations in one end of the 
bay do not effect oscillations in the opposi te'end; the north-south 
basins are "uncoupled." 

Raines [Il967l] analyzed twelve^ selected long-period wave trains 
from three years of tide data (1964-1966) obtained with the Naval Post- 
graduate School (NPS) recorder located on Municipal Wharf #2, Monterey 
harbor. Using graphical methods, he found mean periods of 19-39 min 
and 1. 5-2.0 min. Raines did not rule out the possibility that the 
longer period waves were bay seiching, but he strongly suggested, 
based on his spectra and barometric oscillations occurring simultane- 
ously with half of the wave trains analyzed, that the longer period 
waves were progressive waves produced by air-pressure fluctuations. 

Lacfe of adequate barometric pressure data precluded a correlation of 
the wave periods and air-pressure fluctuations. The shorter period 
waves were judged to be either harbor seiching or surf beat. This 
conclusion was based on the similarity of his results with those of 

r 

Wilson L1965T1 and Hopper Cl967^ . 
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II. INSTRUMENTATION 



The Monterey data for this study were obtained with a standard Coast 
and Geodetic Survey automatic tide gage i Manual of Tide Observations, 
1965H . The gage, which is maintained daily by NPS personnel, is lo- 
cated on Municipal Wharf #2, Monterey harbor. 

This instrument senses changes in water level by means of a float/ 
pulley arrangement. The recording drum is advanced by a clock mecha- 
nism. The drum speed is designed to be 1 in/hr but the NPS gage has 
an hourly feed of approximately 1.06 in/hr. The marigram is recorded 
in rectilinear coordinates on plain white paper. The majority of 
higher frequency wind waves (periods of 4 min and below) are filtered 
by a stilling well which is a 12-in diameter steel pipe with a 1-in 
orfice in the bottom. 

The Santa Cruz data were obtained with a Bristol Model 28 gas- 
purging pressure (bubbler), portable tide'gage located on the Santa 
Cruz wharf. This instrument senses changes in water level by means 
of a nitrogen-filled tube which is connected to a bellows system 
CManual of Tide Observations, 1960 . The marigram is recorded on 
a mechanical -clock strip- chart recorder in curvilinear coordinates. 

The design chart drive speed is 1 in/hr, but a substitution of 
drive gears increased this to 6 in/hr 

A bubbler orfice chamber was connected to the end of the sensing 
tube to reduce wave action. There was also a bellows inlet needle 
valve which could be throttled to further filter wave action from 
the record. The combination of these two filtering mechanisms 
proved inadequate. A stilling well was designed and installed to 
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attain the desired filtering of high frequency wind waves. The well was 
constructed of a 20- ft section of polyvinyl Chloride (pvc) 6-in inside- 
diameter pipe. The well was capped and 16, 1/4-in inside-diameter holes 
drilled in the side. Copper sleeves were inserted to eliminate fouling. 
The 16 holes provided the capability of increasing the orfice from a 
1/4-in to a 1-in diameter opening. 

The response characteristics of the well were determined theo- 
retically for two orfice sizes and three different wave frequencies 
using the equation for the rate of water rise in a well CDoodson and 
Warburg, 194Q : 

rate of rise of water in well = dh.j/dt = 0.6 a/AV2g(ho - h-j ) 
where 

a = orfice area 
A = well area 

0.6= empirical orfice flow coefficient 
g = acceleration due to gravity 
ho = water height outside the well, i.e. 

the forcing function 
h^* = water height inside the well. 

The forcing function, h^, was chosen to be a simple sine function of 
unit amplitude and frequency equal to the wave frequency of interest. 

The initial conditions were h.j = o at t=0. The results are sum- 
marized in Table VI. 

The response characteristics for the Monterey stilling well were 
not calculated since the orfice area to well area ratio is larger, 
providing response characteristics better than those of the Santa 
Cruz well. A 1/4-in orfice was used in the Santa Cruz well. 
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PERIOD 


ORFICE DIAMETER 
(in) 


PHASE LAG 
(deg) 


AMPLITUDE REDUCTION 
{%) 


20 sec 


0.25 


180 


95 


20 sec 


1.00 


72 


45 


60 sec 


0.25 


75 


91 


60 sec 


1 .00 


30 


5 


25 min 


0.25 


6 


1 


25 min 


1.00 


0 


0 



TABLE VI 



Computed Response Characteristics 
of Santa Cruz Stilling Well 
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III. DATA ANALYSIS 



A. SELECTION OF DATA 

Because this study was concerned with the character of seiching 
motions affecting the entire bay (and not their frequency of occurrence), 
only those records with similar long-period wave characteristics at each 
station for equal intervals were considered for analysis. Also, the 
seiching had to persist long enough to provide an adequate number of 
data points for a meaningful analysis. The minimum persistence time 
considered was ten hours. Two days were selected which met these 
criteria; 23 January and 20 April, 1969. 

The Santa Cruz marigram for 23 January was recorded at a rate of 
1 in/hr with a 20- ft range instrument. Although the stilling well 
had not been installed, the trace was relatively narrow and free of 
wind-wave noise. The long-period waves lasted for approximately 15 
hours (1000 23 January - 0100 24 January). 

The Santa Cruz marigram for 20 April was recorded at a rate of 
6 in/hr with a 10-ft range instrument. The stilling well had been 
installed by this date. (In February, the 10-ft range instrument was 
substituted to improve the resolution of the record and the recording 
rate was increased to facilitate digitizing procedures.) The long- 
period waves for this day lasted approximately 10 hours (1100-2100 
20 Apri 1 ) . 

B. DIGITIZING PROCEDURE FOR MONTEREY DATA 

The Monterey record, a rectilinear trace, was digitized using a 
Calma Co. Model 480 mechanical digitizer (Appendix A). The sampling 
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rate, At, for the 23 January record was 33.87 sec whereas the sampling 
rate for the 20 April record was 33.77 sec. The difference in sampling 
rates was not due to the digitizer, which had a constant sampling in- 
terval of O.Ol in, but rather was a result of the different recording 
rates of the two records. The average recording rates were 1.063 in/hr 
and 1.066 in/hr for the 23 January and 20 April records respecti vely. 

C. DIGITIZING PROCEDURE FOR SANTA CRUZ DATA 

The Santa Cruz records were not digitized by the machine method 
because of the problems encountered in digitizing a curvilinear record 
with a rectilinear device. A simple geometrical relationship can be 
used to convert the rectilinear digitized data of a curvilinear 
record, but this conversion introduces errors and high frequency noise 
QSteele, 1967Zlwhich were thought to be excessive for this study. 

These records were, therefore, digitized by hand. It was not possible 
to equal the Monterey data sampling interval of 0.01 in for the 23 
January Santa Cruz record since the latter was recorded at 1 in/hr. 

A sampling interval of approximately 0.025 in was used, giving a 
sampling rate of 89.64 sec. This meant that the sampling rate for 
the Monterey record was about two and one-half times better than the 
Santa Cruz record for this date. 

The sampling interval for the 20 April record, a 6-in/hr record- 
ing rate, was selected to be as close as possible to that for the 
Monterey record. The resulting sampling rate was 33.14 sec, approxi- 
mately one-half second slower than that for the Monterey record. 

The problem, and avoidance, of aliasing in spectral analysis of 
a finite, discrete, record is thoroughly discussed in the literature. 

It is only important to note here that the sampling rates achieved in 
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this study are more than adequate since both tide gages are long-period 
recorders designed to filter out all waves with periods below 4 min. 
Aliasing, therefore, is not considered a problem and is not discussed 
below. 

D. ANALYSIS 

Before the data were analyzed, it was necessary to remove the 
trend due to tidal components. The detrending was accomplished to 
an adequate degree by forming a pure cosine curve, , which closely 
approximated the tidal record for the specific interval digitized, 
and subtracting this "tidal curve" from the raw data, y^- . That is, 

= Acos(2 77fi At-0) i = l ,2,3,. . . ,N 

where 

A =amplitude measured on marigram 
f =tidal frequency measured on marigram 
0 =appropriate phase lag required to match 
data and 7^^ 

N =number of data points, 
and 

Yi = Yi - (m +7|.) i = l,2,3,...,N 

where m is the mean measured on the marigram and Y is the detrended 
data. 

The Fourier coefficients of the detrended data were calculated 
using the IBM/360 library subroutine RHARM. This subroutine is a 
one-dimensional fast Fourier transform (FFT) analysis based on the 
Cooley and Tukey0965l] algorithm. The program is designed to 
analyze N data points where 

N = 2"^ m = 3,4,5,... ,20. 
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The FFT not only greatly reduces the number of calculations from 
earlier Fourier analysis schemes, but also reduces the round-off errors 
in the coefficientSo Specifically, both the number of computations 
and round-off errors are reduced by essentially log2(N)/N [[Cochran, et 
al , 1967Z1 . Before the spectrum was formed, each Fourier coefficient 
was banned. 

The raw Fourier transform is exact for the specific frequencies 
of the calculated coefficients. Adding, however, a term of the form 
Dcos(ft) where the frequency f is not one of the discrete frequencies, 
f j , will alter all the Fourier coefficients in the time series. This 
is, in effect, an energy leakage into the discrete frequencies of the 
raw periodogram. Leakage can be defined as the altering of a spectral 
estimate for a specific bandwidth by the energy at more or less random 
frequencies outside the bandwidth of interest. The effects of leakage 
in the raw periodogram decay as 1/ |f-f-j| as f^ recedes from f. 

Hanning the coefficients before forming the spectrum increases the de- 
cay to l/|f-fil^ [Bingham, et al , 1967^. If a^- and b^. are the raw 
Fourier coefficients, the banned or modified coefficients are formed 
using: 

Ak = -0/4)ak_i + 0/2)ak - (l/4)ak+i 
Bk = -(l/4)b|,_^ + (l/2)b|^ - (l/4)bk+i 
The spectral estimates, (A^ + k=l ,2,3,.. , ,N/2, were normalized 

by multiplying each by the time interval analyzed. 

Figures 17-22 (Appendix B) are the computed spectra for the two 
days analyzed. The 15-hr period wave activity on 23 January is 
depicted by two, overlapping spectra for each station. This is a 
consequence of the RHARM restrictions on the number of data points. 
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The spectra are not the entire spectra, but rather only the low 
frequency portion. The dashed line on the low end of the spectrum in- 
dicates that all energy in the tidal components has not been removed. 
Similarly, the dashed line at the higher end of the spectrum designates 
the start of the high-frequency noise. There is approximately three 
hours difference in the Monterey and Santa Cruz time intervals for the 
spectra for 23 January, which is a consequence of the number of data 
points and the unequal sampling rates. 

The spectra for different stations on the same day show a certain 
amount of correspondence with each other, but it is difficult to match 
periods of oscillation between spectra since it is not easy to deter- 
mine with confidence which peaks in the spectra are spurious and which 
are real . 

Spectral estimates of a random process have a chi-squared distri- 
bution [Bartlett, 1955[ and are inconsistent estimates of the power 
spectrum. Each coefficient has 2 degrees of freedom, and the con- 
fidence in each spectral estimate can be increased by (1) averaging 
spectral estimates over a span of frequencies or (2) "blocking" the 
record, analyzing each block and averaging spectral estimates at 
equal frequencies [Jones, 1965; Hinich and Clay, 1968D. 

The spectra in this study were averaged over three bandwidths. 
Smoothing the spectra in this manner gives more confidence in the 
true peaks, but the resolution in each spectrum is simultaneously 
degraded. The two spectra for each station on 23 January have been 
averaged together at equal frequencies to give a spectrum for the 
entire 15- hr period. 

When sequentially averaging over three bandwidths, there are 
three possible starting points, that is, three possible methods. 
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The spectra were averaged for all three possibilities which permitted de- 
termination of the resolution and stability of each peak. Schematic 
representations of the smoothed spectra are presented in Figures 2-5. 
These spectra represent a compilation of the three averaging methods 
for each spectrum. Only those peaks which were stable in shape and 
energy density are labeled with the bandwidth resolution. All peaks 
not labeled will not be considered further. 

Spectra for 10-hr periods both before and after the long-period 
wave activity on 23 January and 20 April were also calculated for each 
station (Figures 23-30, Appendix B). These spectra were smoothed in 
the same manner described above. Schematics of the "before" and 
"after" spectra for each station are presented in Figures 6-9. 

The original intention in this study was to compute the cross- 
spectra for the Monterey-Santa Cruz records. It was hoped that, with 
the cross-spectral values and average phase lags for each discrete 
frequency, a thorough understanding of the seiching motions in 
Monterey Bay could be obtained. The cross-spectra calculations were 
not feasible, however, for two reasons. 

The records were not accurately synchronized. It is estimated 
that an initial time difference of as much as 3 min between the two 
records could be present. This was due to the variability in the 
seperate gage clocks and the maintenance of the gages by different 
persons at varying intervals. This factor alone is not prohibitive 
since Grauzinis D968I] has outlined a procedure for computing cross- 
spectra for records without initial synchronization. A more serious 
error was created by the unequal digitizing rates between records 
from the two stations. A cumulative phase error was introduced by 
sampling at unequal rates. The data obtained in this study, there- 
fore, are considered unsuitable for cross-spectral analysis. 
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Figure 3 

Spectrum for Monterey - 1000 23 January to 0100 24 January 
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Figure 4 

Smoothed Spectrum for Santa Cruz -1100 20 April to 2058 20 April 
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Figure 5 

Smoothed Spectrum for Monterey - 1100 20 April to 2048 20 April 
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Figure 6 

Smoothed Santa Cruz Before/After Spectra - 23 January 
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Smoothed Monterey Before/After Spectra - 23 January 
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Figure 8 

Smoothed Santa Cruz Before/After Spectra - 20 April 
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Figure 9 

Smoothed Monterey Before/After Spectra - 20 April 
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IV. DISCUSSION 



A. WILSON'S THREE-DIMENSIONAL MODEL 

Wilson's numerical calculations for the three-dimensional oscil- 
lating characteristics of Monterey Bay are depicted in Figures 10-13. 
Figure 10(a) illustrates the grid points used for the solution. The 
numbers adjacent to each grid point are the water depths in feet. 

The dashed line connecting Point Santa Cruz with Point Pinos is the 
assumed boundary nodal line for the solution. 

The successive figures depict decreasing periods (increasingly 
complex modes) of oscillation. The contours are water level ampli- 
tudes normalized to the highest anti -node for the mode. The inset 
to the left of each figure is a simplified modal oscillation. The 
sequence of periods is 

T_ = 44.2, 29.6, 28.2, 23.3, 21.6, 20.4, 

19.4, 18.7, 17,6 13.3,..., 

12.4,. ..min. 

Wilson notes that the assumed boundary nodal position could 
very well be incorrect and might be more to seaward, closer to 
the 100-fathom curve. Changing the assumed nodal line would alter 
the boundary conditions and would affect both the modal periods and 
the geometry of the oscillations. 

Modes 1 and 2 indicate strong oscillations in the northern 
portion of the bay and little effect in the southern portion of the 
bay. Mode 3 shows the first strong oscillation in the southern 
portion. Mode 4 indicates similar oscillations at both the northern 
and southern ends with weak oscillations over the canyon. Successive 
modes become more complex. 
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Numerical Calculations of Modes of Oscillation of 
Monterey Bay (from Wilson, 1965) 
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Numerical Calculations of Modes of Oscillation of 
Monterey Bay (from Wilson, 1965) 
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Numerical Calculations of Modes of Oscillation of 
Monterey Bay (from Wilson, 1965) 
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Numerical Calculations of Modes of Oscillation of 
Monterey Bay (from Wilson, 1965) 



Wilson does not state the highest mode in which he has confidence, 
but he does note that errors will increase in any eigenvalue problem as 
the mode number increases. This increasing error can be slowed by 
choosing a finer grid system. 

Wilson argued that his observations (Tables IV and V) tend to con- 
firm his computations. Note that for modes exhibiting strong oscil- 
lations in only one end of the bay, amplitudes in the oscillating end 
are 10 to 100 times larger than amplitudes at the opposite end. (Only 
modes 7,8, and 20 show nearly equal amplitudes at the bay extremities.) 
This led Wilson to conclude that oscillations in one end of the bay 
have no effect on, and occur independently of, oscillations in the 
opposite end. 

B. INTERPRETATION OF SPECTRA 

There are several obvious statements which can be made about the 
smoothed "before", "during", and "after" spectra. For ease in com- 
paring the energy levels of the various spectra, all the spectra for 
the same day at both stations are plotted to the same scale in Figures 
14 and 15, These schematics show only the smoothed energy peaks. 

The energy levels for the during spectra are always higher at 
Santa Cruz than at Monterey, but the energy levels for the before/ 
after spectra for both stations on both days are the same. The 
energy levels for the during spectra at Santa Cruz are much greater 
(23 January) and greater (20 April) than the energy levels in the 
before/after spectra for Santa Cruz, In Monterey, however, the 
energy levels in the during spectra are only slightly more (23 
January) and slightly less (20 April) than the energy levels in the 
before/after specta. No convincing explanation can be offered for 
these differences. 
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FIGURE 14 

Schematic Representation of Spectra for 23 January 
Santa Cruz and Monterey 
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Schematic Representation of Spectra for 20 April 
Santa Cruz and Monterey 
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There is evidence in the during spectra that 16-18, 13-15, and 45-60 
min oscillations were present at both stations for the two days analyzed. 
About 10 other spectra for each station were calculated in the course of 
this study. The 16-18 min oscillations were consistently present where- 
as the 13-15 min oscillations were usually, but not always, present. 

The before/after spectra all exhibit some type of long-period wave 
activity. Most show oscillations in the 45-60 min range. There does 
not seem to be any correspondence between the before/after spectra for 
different stations on the same day. It is possible, for some periods, 
to detect signs of build-up and decay in the before/after spectra at 
a station. On the other hand, there are some periods which are present 
in only the before or only in the after spectrum and not present in the 
during spectrum. 

Before attempting to analyze the spectra, it is well to ask if 
the long-period waves on 23 January and 20 April were seiching in the 
bay. Unfortunately, it is impossible to tell by simply examining the 
power spectra. An accurate determination could be made if the records 
had been accurately synchronized so that phase angles for each fre- 
quency could be calculated or if data from an off-shore recorder were 
available so that the off-shore and near-shore spectra could be com- 
pared. Without this information, however, one can only say that the 
two most distinct possibilities for waves in this frequency range 
are (1) bay seiching or (2) shelf oscillations in the vicinity of 
the recording instruments. 

Wilson's model for seiching in Monterey Bay predicts greatly dif- 
ferent amplitudes at the extremities of the bay, but of course the 
periods of oscillation are the same over the entire bay. MunklZl962ZI 
found shelf oscillations ranging in period from 2.5 min to 3.1 hr. 
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These osci 1 Idtions are aetermined by the dimensions of the continental 
shelf and extend from the shelf edge to the shoreline. The oscillations 
might be quite localized with different periods of oscillation for rela- 
tively near points along a coast. Munk notes that shelf oscillations can 
radically shape the spectrum of a record measured near shore, since the 
anti-node is at the shoreline. It is possible that bay seiching and 
shelf oscillations can occur at the same time. The record from an off- 
shore recorder would be invaluable in defining the periods of shelf 
oscillations. 

Notice for the 23 January spectra that Santa Cruz has a strong 
oscillation at 27.0-27.6 min and Monterey has a strong oscillation at 
32.1-38.5 min. These two oscillations are not repeated on the 20 April 
spectra. It could be argued that these are localized shelf oscillations 
as found by Munk. It could also be argued, however, that these are 
periods of bay seiching and that the amplitudes of these periods were 
so small on the opposite end of the bay that they did not show in the 
spectral analysis, as predicted by Wilson's model. Several other 
periods of oscillation in the spectra can be treated in the same manner 
as this, but they are not as pronounced as these two. However, it is 
not possible to take every peak in the spectra for 23 January and 20 
April and state that they are evidence of bay seiching. 

It is more reasonable to look for periods of oscillation which are 

evident at both stations on both days analyzed. If these occur and 

have similar amplitude ratios on the two seperate days, then they are 

probably periods of bay seiching. Two bandwidths were found: 15.7-21.5 

min (hereafter referred to as the average, 18.6 min) and 12.8-14.1 min 

(13.4 min). All the during spectra have oscillations within these two 
> 

bandwidths. These two periods then are selected as being associated 
with bay seiching. 
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Define an amplitude coefficient, A, to be: 



A = amplitude at Santa Cruz 
amplitude at Monterey 

The amplitude coefficients for the 13.4 and 18.6 min oscillations can 
be easily computed. Multiplying the energy density in a bandwidth by 
the bandwidth will give a measure of the energy. Taking the square 
root of this energy gives a relative amplitude for each station. The 
amplitude coefficients for the 13.4 and 18.6 min oscillations are 
given in Table VII. The coefficients are in good agreement, and are 
much less than would be predicted by Wilson's model. 



TABLE VII 

Amplitude Coefficients for 23 January and 20 April 



BANDWIDTH 


AMPLITUDE 


COEFFICIENT 


(min) 


23 January 


20 April 1 


15.7-21.5 


2.78 


1.74 


12.8-14.1 


2.20 


3.80 



Thus, the two periods of oscillation selected in this report 
have similar characteristics and probably represent seiching in 
the bay. Also, since the amplitudes at Monterey and Santa Cruz 
are of the same order of magnitude, it appears that oscillations 
can occur which do not act independently of each other with great- 
ly different amplitudes at each end. This is not to say that the 
canyon cannot be a nodal line for bay seiching (it is certainly 
reasonable that oscillations over the canyon are quite small, as 
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predicted in Wilson's model), but it is not reasonable, based on the 
above results, that the canyon uncouples all seiching in the northern 
and southern halves of the bay. 
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V. SUMMARY 



The Monterey Submarine Canyon has an effect on the seiching motions 
in Monterey Bay, but does not appear to divide the bay into two, in- 
dependent oscillating basins. Some seiching action in the bay on 
23 January and 20 April 1969 occurred over the entire length of the 
bay, having amplitudes of the same order of magnitude at Monterey and 
Santa Cruz. 
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APPENDIX A 



Use of the Calma Co. Model 480 Digitizer 



The Calma Model 480 digitizer located at the Fleet Numerical 
Weather Central' (FNWC) Point Pinos annex was used to digitize the 
Monterey data.’ This method of digitizing proved to be both faster 
and more accurate than template methods. 

The digitizer reads and records data in the X and Y directions 
with a sampling interval which can be set to 0.01, 0.02, 0.04, 0.08, 
or 0.15 in. The sampling interval used for this study was 0.01 in 
for both the X and Y directions. The maximum absolute sampling error 
for the machine is 0.012 in. The output is external BCD, stored on 
556 bpi , 7 channel, tape. The tape can be made compatible with the 
IBM/360 system. It was easier, however, to use the CDC/6500 computer 
at FNWC. 

A simple program, CONVERT, for interpreting the digitized record 
is presented in Figure 16. This program accomplishes six steps: 

1. Converts the data from external BCD to display code. 

2. Arranges the data in column matrices of 80 character length. 

3. Interrogates each character in the matrix to determine if the 
character is a 

a. flag. 

b. plus or minus travel in the X-direction. 

c. plus or minus travel in the Y-di recti on. 

4. ..Sums the Y-direction travel. 

5. Prints the Y value, V(m), for each 0.01 in increment travel 
in the X- direction. 

6. Punches cards for the amplitude values which were subsequently 
used as data input for RHARM on the IBM/360. 

Headings may be entered on the tape with the keyboard control, 
but CONVERT is not designed to read tape headings nor does the pro- 
gram have the capability to search for a particular set of data on 
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the tape. Consequently, all data on the tape will be analyzed each 
time CONVERT is used. This requires that the user know the sequence 
of digitized data on the tape. 

The digitizing procedure for which CONVERT is designed is listed 
below. This procedure is followed for each interval digitized. 

1. Flag the tape using the "Flag" foot pedal on the machine. 

2. Digitize the interval desired. 

3. Flag the tape twice. 

4. Enter IRG on the tape by pressing the proper key on the keyboard 
console. 

Detailed information on the operation of the digitizer, plus limited 
computer programming information, can be found in the Model 480 
Digitizer Instruction and Maintenance Manual. 
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PROGRAM CONVFRT ( I NPUT, OUTPUT , PUNCH) 



OIMFNSION U( 500-^) ,V(500'^) , N( 80 ) ,NK( 80 ) ,IBUFF(300n ) 

READ 1CC,L 

100 FrRMAT(I2) 

on 125 JJ=1,L 
READ 101,AVGX 

101 FORMAT(FIO.O) 

0CLT=36.C/AVGX 

C AVGX IS THE AVERAGE DISTANCE THE RECORDING DRUM OF THE 

C TIDE GAGE ADVANCES PER HOUR. 

PRINT 1C2,DELT 

n2 FnRMAT(lHC,2X,25H, SAMPLING INTERVAL E QU AL S, I X , E 1 5 . 7, 
11X,7H SECONDS) 

PRINT 103, JJ 

103 FORMAT ( IHO ,9H JJ EQUAL S , 2X , I 2 ) 

DO ICA 1=1,5000 
U( I )=0 .0 
lOA V(I)=O.C 

DO 1C6 1 = 1 ,3000 

106 I BUFF( I )=0 
CC’UNTX=0.'^ 

COUNTY=C.C 

NUM=3C0C 

M=1 

K8=0 

CALL L lOF ( 5LRBCD1 , I BU FF, NUM, NPAR , NEUF ) 

IF(NEOF) 602,107 

107 K=-7 
KF=0 

108 K=K+8 
KA=K+8 
KC=0 

DO 109 KB=K,KA 
KC =K C + 1 

109 NK(KC)=IBUFF(KB) 

OPCODE (80, 110, NK) ( N ( I ) , I =1 ,80 ) 

110 F0RMAT(8CR1) 

DO 12C 1=1,80 

IF(N( I ).E0.47B ) GO TO 111 
IF(N( I ).EQ.50B) GO TO 113 
IF(N( I ).E0.55B) GO TO 114 
IF(N( I ) .EG .34B ) GO TO 115 
SYMBOL ''«”,(47B), IS A FLAG. 

SYMBOL "/", (50B), REPRESENTS AN INCREMENT TRAVEL IN 
THE MINUS X OR Y DIRECTION. 

SYMBOL "C", (55B), REPRESENTS ZERO TRAVEL IN THE X OR 
Y DIRECTION. 

SYMBOL "I”, (34B), REPRESENTS AN INCREMENT TRAVEL IN 
THE POSITIVE X OR Y DIRECTION. 

GO TO 120 

111 PRINT 112, M, I 

112 FORMAT(1HO,2I1C) 

IF(M.GT.IC) GO TO 121 
GO TO 120 

113 RX=-0.01 
K8=K8+1 
GO TO 116 

114 RX=0.0 
K8=K8+1 
GO TO 116 

115 RX=0.01 
K8=K8+1 

116 K3=K8/2 
K3=2*K3 

IF(K3.EQ.K8) GO TO 118 
COUNTX=COUNTX+RX 
IF(COUNTX.NE.O.O) GO TO 117 
GO TO 12C 
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117 U(M)=COUNTX 
GO TO 12C 

118 CCUNTY=CGUNTY^RX 
IF(COUNTX.NE.O.O) GO TO 119 
GO TO 120 

119 V(M)=CCUNTY 
COUNTX=0.0 
M=M+1 

IF(M.GT.5CC0> GO TO 600 

120 CONTINUE 

121 TIME=(M*DELT)/3600. 

PRINT 122, TIME 

122 FORMATdH ,1CX,28H TOTAL TIME OF RECORD EQUALS, 2X, 
E15.7,2X,7H HOURS. ) 

PRINT 123, IV(1 ),I=1,M) 

123 FORMATdH ,10X,14F7.2» 

PUNCH 12A, (V ( I ) , I = 1,M) 

124 FORMAT (14F5. 2) 

125 CONTINUE 
STOP 

600 PRINT 601 

601 FORMAT (IHl ,20X , 37H ***♦ U AND V SPACE INADEQUATE ****) 

602 STOP 

END 
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